Heat release in clusterization of nanoparticles 
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We show that the effect of size dependence of the melting temperature of nanocrystals may be 
used to govern anomalies of thermodynamic properties of nanocrystals in the premelting range. For 
example, if temperature of nanocrystals is near the melting point, their mechanical contact should 
result in a change of the temperature and release of some heat. We discuss fundamental possibility 
of construction of a "heat pump" to extract thermal energy from a low-temperature medium. 

PACS numbers: 64.70.D-, 65.80.+n 



Reduction of crystal size up to the nanoscale regime 
is always accompanied by modification of its phonon 
and electronic spectra, thus changing both electrophysi- 
cal and thermodynamic properties of the system. Such 
effects may be so dramatic, that the state of the matter 
would change fundamentally. 

Nanocrystals display strong size dependence of ther- 
modynamic properties. For example, heat capacityii^i^, 
Debye temperature^, and melting temperature^i^i^ vary 
with particle size even in the mesoscopic size range. As 
shown in Ref. such behaviour of thermal properties of 
nanocrystals is caused by size-dependent rearrangement 
(quantization) of their phonon spectra. Such modifica- 
tion of the vibrational spectrum results in changes of sta- 
tistical parameters of atoms composing the nanocrystal 
and, therefore, in size-dependent changes of average val- 
ues of energetic contributions to the free energy of the 
system. 

A temperature Tc of anharmonic instability of the crys- 
tal is especially sensitive to the above mentioned rear- 
rangement of the eigenfrequency spectrum due to reduc- 
tion of the crystal size^. Approaching is accompanied 
by nonlinear rise of some thermodynamic properties^ 
and dramatic reduction of the formation energy of crys- 
tal structure defect o^°'^^i^^ . These facts allow one to re- 
late Tc with the melting temperature Tm and to con- 
sider anomalous behaviour of thermodynamic properties 
as premelting effects^. For nanocrystals, it is important 
that the instability point can be approached not only by 
increasing temperature, but via reduction of the crystal 
size, too. 

To determine thermodynamic properties of nanocrys- 
tals in Refs.FzlfTsl a self-consistent statistical method® was 
used. The method had been developed to calculate ther- 
modynamic characteristics of simple solids in wide range 
of temperature and pressure using statistical distribution 
functions of atomic displacements. Parameters of these 
distribution functions were determined by phonon spec- 
trum only. Size influence on parameters of such statisti- 
cal distribution for a nanocrystal was taken into account 
by consideration of size-dependent rearrangement of the 
vibrational spectrum''. The nearest- neighbour interac- 
tion was approximated with the Morse potential with 
parameters determined by fitting the observed values of 
interatomic distance, bulk modulus, and sublimation en- 
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FIG. 1: Melting temperature of a free Cu thin plate versus 
thickness (by Ref. 0). 



ergy at T = 

In Fig. [T] we present the melting temperature of a 
Cu nanocrystalline plate versus thickness h computed in 
Ref. '7 from the assumption that Tm ~ Tc- Similar size- 
dependent melting curves are inherent to all kinds of free- 
standing nanocrystals, e.g. rare-gas nanocrystals^'^, thus 
implying universality of the mechanism of size-dependent 
quantization of vibrational spectra which is responsible 
for a shift of temperature dependence of thcrmodynami- 
cal properties of nanocrystals. 

Along with the melting point, premelting anomalies of 
thermodynamic properties of nanocrystals also shift to- 
wards lower-temperature range. This effect is illustrated 
with Fig. [2] which displays temperature dependence of 
isobaric heat capacity Cp of Cu nanoplates of different 
sizes calculated from the equilibrium value of the crystal 
free energy^. A dramatic increase in Cp corresponds to 
the premelting temperature range, in accord with results 
of calorimetric studies of nanocrystalline materials^'^. A 
similar behavior of Cp{h) was also theoretically found for 
the rare-gas nanocrystals^'^. As crystal size decreases, the 
premelting anomaly takes place at lower temperatures. 
Thus, size of a nanocrystal may be considered as an in- 
dependent thermodynamic variable, in addition to pres- 
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FIG. 2: The temperature dependence of the isobaric heat ca- 
pacity of Cu nanocrystals: (a) bulk crystal; (b) h — 10.2 nm; 
(c) h = 5.1 nm; (d) h = 2.55 nm; (e) h = 1.26 nm (by Ref.Q)- 



FIG. 3; A change in the specific heat of fusion of a nanocrys- 
talline plate versus thickness calculated for Cu and Ag. 



sure and temperature. It is important that variation of 
the crystal size allows one to govern the premelting tem- 
perature range, where thermodynamic properties show 
anomalous behaviour. 

The above size-dependent shift of the premelting range 
leads to a number of physical effects that should be ob- 
served in nanocrystals. For example, a change in the 
specific heat of fusion A of a nanocrystal of size h in 
comparison with the bulk value can be approximately 
represented as 



AA = e(r„,(/^),;i)-£"(r^), 



where 
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is internal energy per mass unit, g{T, P) is Gibbs free 
energy per mass unit, superscript denotes bulk values of 
thermodynamic parameters. Using an explicit expression 
for g{T, P) derived in Ref. 0, it is easy to show that AA 
depends on h as 
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where 7 is a factor taking into account discrete character 
of the vibrational spectrum of the nanocrystal, R is inter- 
atomic distance. In Fig. [3] we plotted AA of a nanocrys- 
talline plate as a function of its thickness calculated for 
Cu and Ag. 

We may conclude from Fig. [2] that an excess of spe- 
cific heat accumulated in nanocrystals of smaller size 
should be released at tight mechanical contact of a num- 
ber of particles, which is equivalent to formation of a 
larger particle. Such heat release should lead to a tem- 
perature change of nanocrystals and then to a heat ex- 
change between particles and their surrounding. For 



quantitative description of such an effect, let us consider 
free nanocrystalline plates of thickness h ^ R. If two 
such nanocrystals are brought into tight mechanical con- 
tact, the vibrational modes with wavelengths larger than 
a characteristic thickness of the contact become collec- 
tivized, and their frequencies become equal to the fre- 
quencies of corresponding modes of a crystal with thick- 
ness 2h (see, e. g., Ref. ,15). Breaking the contact, i.e. 
separating the two thin plates, results in the reverse re- 
arrangement of their vibrational spectra and is accom- 
panied by a change Ash in the specific internal energy. 
Obviously, an energy change due to joining the plates 
is — Ae^. If the nanoparticles are isolated, the energy 
change would lead to a change in the temperature of the 
nanocrystals, AT — Ash/cv, where cy is specific iso- 
choric heat capacity. In Fig. [1] we plot a contact temper- 
ature change AT of two thin plates versus reduced ini- 
tial temperature T/T^ calculated for Cu and Ag plates 
of thickness h = 5.1 nm {Tm{h)/T^ = 0.919 for Cu and 
0.914 for Ag). As follows from the results presented in 
Fig. m AT is maximal if the initial system of separated 
nanoparticles is close to the melting point. Fig. [5] rep- 
resents temperature change due to contact of two non- 
metallic (Xe) nanocrystals computed for two values of 
the plate thickness, h = 5.6 nm and h — 3.5 nm. 

A contact temperature difference AT between a 
nanocrystalline cluster and surrounding medium should 
initiate transport of some heat from the nanoparticles to 
the medium. We can estimate a transferred amount of 
heat from the following reasoning. Heat amount accumu- 
lated by a particle of mass m and size h at temperature 
T is 



Q — m cp(t, h) At, 
Jo 



(4) 



where cp{T, h) is isobaric specific heat. Usually, a change 
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FIG. 4: A temperature change at mechanical contact of two 
metalhc thin plates of thickness h = 5.1 nm versus initial 
temperature calculated for Cu (theoretical value of the bulk 
melting temperature = 1425 K) and Ag (T^ = 1228 K). 



FIG. 6: Amount of transferred heat versus initial temper- 
ature calculated for Cu and Ag nanocrystals of thickness 
h — 5.1 nm. 
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FIG. 5: A temperature change at mechanical contact of two 
Xe thin plates of thickness h versus initial temperature cal- 
culated for two values of h. 

in Q is ascribed to a change in temperature of the 
medium. However, in the case of a contact of nanocrys- 
tals whose initial temperature is close to the melting 
point Tra, there is a remarkable change of cp (Fig. [2]), 
so that 

AQKim [ [cp{T,2h) - cp{T,h)] AT. (5) 

Obviously, the value of AQ increases substantially as the 
nanocrystals approach the melting point (Figs. [H [7]). For 
example, an estimate of the heat released in the contact of 
Cu plates of thickness h = 5.1 nm near Tm{h) — 1310 K 
gives about 12 J/g. 



If AT - 10 - 100 K, the effect of contact change of 
temperature can be used for mechanical heat transport 
from media with a lower temperature to media with a 
higher temperature. Indeed, transfer of nanocrystals of 
size h from medium 1 with temperature Ti into medium 2 
with temperature T2 > Ti and their clusterization in the 
medium 2 results in heating of the nanocrystals and fur- 
ther heat transport from the nanocrystals to the medium 
2. On the next step, returning the clustered nanocrys- 
tals to the medium 1 and their separation leads to cooling 
of separated nanocrystals and, therefore, to heat trans- 
port from the medium 1 to the nanocrystals. Realization 
of such a cyclic process would allow one to construct a 
"heat pump" to extract thermal energy from lower tem- 
perature media. It is important that, due to mechanical 
work performed over the system, this scheme does not 
contradict to the second law of thermodynamics: "Heat 
generally cannot flow spontaneously from a material at 
lower temperature to a material at higher temperature" 
(Clausius statement). 

The effect of heat release and change of temperature 
due to contact of nanoparticles may play important role 
in biological processes, e.g. as a channel for viruses to 
obtain energy from the environment. 

To conclude, it is important to emphasize that the 
present study is focused on consideration of thermal ef- 
fects in nanocrystals of mesoscopic size range, so we can 
neglect effects due to structural changes inherent in nan- 
oclusters. Such structural transitions are observed, e.g., 
in Au nanoclusterai^iilii^ and result in hysteresis phe- 
nomena at melting. It should be also noted that size- 
dependent thermodynamic effects in nanocrystals are 
caused by quantization of phonon spectra and have uni- 
versal nature. Contribution of conduction electrons to 
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size effects is proportional to a factor of Xp/hr^, where 
Ai? is electron wavelength at the Fermi level. For met- 
als, Xp/h <^ 1, therefore, electronic influence on size- 
dependent changes of thermodynamic properties can be 
neglected. This contribution may be important, e.g., 
in nanocrystallinc degenerated semiconductors^^ where 
Xp/h < 1. 
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FIG. 7: Amount of transferred heat versus initial temperature 
of Xe nanocrystals calculated for two values of h. 
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